Random— phase approximation for the grand— canonical potential of composite 
fermions in the half— filled lowest Landau level 



00 

PL, 
(N 



J. DieteP, Th. Koschny Q , W. Apel & , and W. Weller a 
a Institut fur Theoretische Physik, Universitat Leipzig, 
Augustusplatz 10, D 04109 Leipzig, Germany 
b Physikalisch-Technische Bundesanstalt, 
Bundesallee 100, D 38116 Braunschweig, Germany 
(February 1, 2008) 

We reconsider the theory of the half-filled lowest Landau level using the Chern-Simons formu- 
lation and study the grand-canonical potential in the random-phase approximation (RPA). Cal- 
culating the unperturbed response functions for current- and charge-density exactly, without any 
expansion with respect to frequency or wave vector, we find that the integral for the ground-state en- 
ergy converges rapidly (algebraically) at large wave vectors k, but exhibits a logarithmic divergence 
at small k. This divergence originates in the k~ 2 singularity of the Chern-Simons interaction and it 
is already present in lowest-order perturbation theory. A similar divergence appears in the chemical 
potential. Beyond the RPA, we identify diagrams for the grand-canonical potential (ladder-type, 
maximally crossed, or a combination of both) which diverge with powers of the logarithm. We 
expand our result for the RPA ground-state energy in the strength of the Coulomb interaction. The 
linear term is finite and its value compares well with numerical simulations of interacting electrons 
in the lowest Landau level. 
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I. INTRODUCTION 

The combination of an electronic interaction and a 
strong magnetic field in a two-dimensional electron sys- 
tem yields a rich variety of phases. These are best clas- 
sified by the filling factor v, which is the electron den- 
sity divided by the density of a completely filled Landau 
level. In the case of v = 1/2, the behavior of the sys- 
tem resembles that of a Fermi liquid in the absence of a 
magnetic field, or at small magnetic fields, and over the 
past years, a intriguing picture has emerged : at v = 1/2, 
each electron combines with two flux quanta of the mag- 
netic field to form a composite fermion; these composite 
fermions then move in an effective magnetic field which 
is zero on the average. The interpretation of many ex- 
periments supports this picture. Here, we mention trans- 
port experiments with quantum (anti-) dotaj, in which 
features of the resistivity are related to closed loops of 
the composite fermions around the dots, and also focus- 
ing expcximentsa; more references to recent work can be 
found irJ3i3. 

Theoretically, the picture was developed in the pioneer- 
ing work of Halperin, Lee, and ReadQ. They formulated 
the Hamiltonian in terms of Chern-Simons (CS) trans- 
formed electrons and studied within the random-phase 
approximation (RPA) many physical quantities. Most 
prominent among these is the effective mass of the com- 
posite fermions which they found to diverge at the Fermi 
surface. This problem of the effective mass was again 
taken up by Stern and HalperirJS who calculated the ef- 
fective mass and also the quasiparticle interaction func- 
tion of Landau's Fermi-liqwid theory, see alscraiJ. Re- 
cently, Kopietz and Castillaliil used the method of higher- 



dimensional bosonization to extend the analysis beyond 
the RPA. Thus treating the infrared singularities, they 
found that the vertex corrections render the quasiparticle 
mass finite. But the composite fermion spectrum should 
scale with the strength of the_Coulomb repulsion. Most 
recently, Shankar and MurthyEj proposed, within the CS 
approach, a new parametrization which allows to sepa- 
rate the contributions of the CS particles from those of 
the magneto-plasmon oscillators in the Hamiltonian. Af- 
ter they restrict the number of the magneto-plasmon os- 
cillators to the number of electrons in the original Hamil- 
tonian, they get a finite quasiparticle mass which scales 
with the inverse of the strength of the Coulomb repulsion. 

The CS formulation of interacting electrons, with a 
density such that the filling factor of a Landau level 
is v = 1/2, nicely realizes the concept of composite 
fermions, in accordance with the experimental observa- 
tions. The price to be paid for using the CS transforma- 
tion is that the non-interacting system of electrons in a 
strong magnetic field, after the transformation, becomes 
a higly complicated, interacting problem - formulated 
in terms of the CS field. The CS interaction diverges 
for small wave vector k as 1/k 2 . Therefore, the singular 
diagrams in leading order need to be resummed; these 
diagramms are well known as the RPA of the jellium 
problem in d — 3. But there is no small parameter in the 
problem of non-interacting electrons in a strong mag- 
netic field and there is no reason to expect that the RPA 
diagrams already reproduce the correct result. In fact, 
as we shall see, the RPA is even plagued with a marginal 
divergence in the energy. The understanding, and even 
the evaluation of the RPA, still seem to be incomplete: 

Is it necessary to introduce an ultraviolet cut off in 
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the momentum integrals, in order to restrict the number 
of the magneto-plasmon oscillators, or is that already 
inherent in the theory? After all, the introduction of 
the cut off seems to be crucial for the conclusions in the 
work of Shankar and Murthy. The simplest quantity for 
a study of this question is the ground-state energy. Its 
evaluation in RPA should be straight forward. But all 
approaches, of which we are aware, concentrate on the 
behavior at small wave vector. Therefore, we reconsider 
the RPA and calculate the expression for the energy, re- 
fraining from any approximations in the bare response 
functions. 

There is a second motivation for studying the energy 
in RPA. For physical reasons, it is clear that the exact 
energy per particle should be_cxpandable in the strength 
of the Coulomb interactionEj. The lowest-order term 
is exactly known, of course, it is the energy of the un- 
perturbed lowest Landau level; still, one can not expect 
to obtain that result (or the corresponding degeneracy) 
correctly within the RPA. Given the lowest order, we can 
proceed calculating - in RPA - the next, the term lin- 
ear in the strength of the Coulomb interaction. Within a 
restriction to wave functions of the lowest Landau level, 
there are no corrections beyond this first order. The con- 
tribution of the Coulomb interaction to the total energy 
seems not to be measurable; still, there are numerical 
simulations of interacting electroaSjin the lowest-Juandau 
level by Morf and d'Ambrumeniji-3 and Girlichta which 
are well suited for a comparison with our analytical re- 
sult. We want to add that such an expansion with respect 
to the Coulomb interaction should be useful also for the 
calculation of other quantities, as the self energy of the 
composite fermions, because there the scale is set by the 
Coulomb interaction. 

The main problem of the RPA, which does not exist 
for interacting electrons at zero magnetic field in d = 3, 
remains the logarithmic infrared divergencies caused by 
the 1/fc 2 behavior of the CS interaction at small wave 
vector. Such a divergence already appears in the lowest 
order in the RPA. Beyond the RPA, there are diagrams 
showing divergencies with powers of the logarithm. 

The outline of this work is as follows. First, we intro- 
duce the model in its CS formulation. Then, we calcu- 
late ground-state energy and chemical potential from the 
grand-canonical potential in RPA using the exact bare 
response functions. Further, we calculate the first term 
of the expansion of the energy in the electron-electron 
(Coulomb) interaction. Finally, we discuss the logarith- 
mically divergent diagrams beyond the RPA. 



II. CHERN-SIMONS DESCRIPTION 

We consider interacting (spin aligned) electrons mov- 
ing in two dimensions in a strong magnetic field B in 
the negative z-direction transverse to the system. The 
electronic density of the system is chosen such that the 



lowest Landau level of a non-interacting system is filled 
to a fraction of v — where 4> is an even integer. We 
are mainly interested in.-j/> = 2. Following the work of 
Halperin, Lee, and React], we describe the Hamiltonian 
in terms of Chern-Simons (CS) fermions or composite 
fermions which consist of the original electrons and cf> 
attached flux quanta : 
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Here, creates (and a k annihilates) a CS fermion with 
wave vector k. Since we consider only cases where the 
filling factor is given by 1/0, the sum of the external 
magnetic fielcLand the mean CS field, AB, vanishes 
(compare Ref.Q) and is given by the free spectrum 
(-k = k 2 /(2m,b), where mj, is the electron band mass; we 
use h = 1, but restore h in final formulae. The chemi- 
cal potential is denoted by fi. Because of AB = 0, we 
can use periodic boundary conditions; F is the area of 
the sample. The interaction has contributions from the 
fluctuations of the CS vector potential and from the orig- 
inal Coulomb interaction of the electrons, V(r) = e 2 /r, 
where e 2 = e 2 / (Ane r eo) an d £ r is the dielectric constant. 
Collecting all these contributions, we get 

W k;kl)k2 = W(k) + ^ i f(ki - k 2 ) x e z ; (2) 

e 2 denotes the unit vector in z direction, orthogonal to 
the system. Here, the last term comes from the term 
linear in the CS vector potential. The quadratic term is 
contained in K(k) in 



W(k) = K{k) + V(k) , (3) 
K(k) = 

V(k) = 



K(k) + V(k) , 
(2tt0) 2 N 1 

rrib F k 2 ' 
2^e 2 



Rewriting the original Hamiltonian in terms of CS oper- 
ators a) , a as in the Hamiltonian ([j]), we have made an 
approximation which is already implicit in the treatment 
of Ref.Q : In the term quadratic in the CS vector po- 
tential, we have replaced a pair of operators, a^a, by its 
average and neglected thus the third order fluctuations 
in the density, which would lead to three particle inter- 
actions, i.e., sixth order in a', a. This average leads to 
the appearance of the density N/F in the CS interaction 
K(k); N is the mean electron number, the CS particle 
number. For the filling factors v which we consider here, 
we have 
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Ib is the magnetic length, l^ 2 = \eB\/K. The Hamilto- 
nian ([!]) is the starting point for a standard perturbation 
theoretical (diagrammatic) treatment. In the following, 
we will limit ourselves to a study of the grand-canonical 
potential within the RPA. This will enable us to derive 
explicit results for the ground-state energy and the chem- 
ical potential. 



III. GRAND-CANONICAL POTENTIAL 



the RPA without additional approximations. Therefore, 
we need to calculate the unperturbed response func- 
tions exactly. They are given by a one-loop integral 
involving two unperturbed single-particle Greens func- 
tions. The density-density response function (polariza- 
tion part), Iloo, was-jCalculated for the jellium problem in 
d = 3 by LindharcO. In the present case of d = 2, fortu- 
nately, both response functions can again be calculated 
analytically. After performing the sum over Matsubara 
frequencies (of Fermi type), we get 



The calculation of the grand-canonical potential 
Q(T, F, n) in the RPA is described in text booksM The 
only unusual feature in the present case is the (linear) de- 
pendence of the interaction Wk ; ki,k 2 on ki k 2 . Thus, 
one has to consider two kinds of vertices in the pertur- 
bation theory, a density vertex and a current vertex, and 
the RPA equation becomes a matrix equation as stressed 
already in Ref.Q. The grand-canonical potential in the 
RPA, f^RPA, is given by the series of ring diagrams: 




(3n$, A (T,F,ri=0 



(•5) 



Here, we indicate the order of the leading singularity as 
the infrared wave vector cutoff, rj cx decreases. 
For n — 1, this is a logarithmic singularity. The singu- 
larity orginates from the divergence of the CS interaction 
K(k) at small wave vector. The symmetry factor, l/(2n), 
is explicitely written in the diagram. We use the Mat- 
subara technique and write the result as 

n RPA (T,F,fi) = 
-l£ln{l + e^->} 



-4 E E ln i 1 - n °°( fc ' w ) + n n( fc ' ^ 



(6) 



Here, the first term is the grand-canonical potential of 
the free CS fermions. The second term is the contribu- 
tion of the magneto-plasmon oscillators. In (||), the sum 
is on Matsubara frequencies of Bose type (u> = 2wn/(3, 
n integer) and (3 = l/(fceT). The response functions of 
the unperturbed system, n o and IIn, are defined and 
studied in the next section. 



A. Unperturbed response functions 

As explained in the introduction, our aim is to eval- 
uate ground-state energy and chemical potential within 



n 00 (fc,w) = 



d 2 p e- M ° n F {e p -j£) 
(2tt) 2 -iio + e p - £| p _ k | 
(to — ► — to) 



(7) 



and 



n llM = (^) 

m b 



d 2 p (p x k) 2 



(27r) 2 k 4 

n F (e p - n) 



— lU) - 



|p-k| 



(8) 



Here, np{x) is the Fermi function. The functions noo 
and nn depend only on the absolute value of the wave 
vector k. Note that in order to simplify our formulae, we 
absorbed a prefactor in rin. In the limit of zero temper- 
ature, when tif becomes the step function, it is a tedious, 
but straight forward task to perform the remaining inte- 
grals exactly without any approximative expansion in k 
or lu. The zero-temperature results are 



n 00 (fc,cj) = - — 



47T 



,io k . , 
e k 2k F 



-co) (9) 



and 



En (ft, a;) 



TT(. 



48m 6 

+ {U1 -> -LO) , 



(10) 



where kp is the Fermi wave vector defined by /i 
k 2 F /(2rrib)- We use the abbreviation 



z(x, u) = 1 + ix — v(l + ix) 2 



(11) 



and we define the square root as analytical continuation 
from positive values of the root for positive radicands. 
Our results, the zero-temperature response functions noo 
and nn, are correct for arbitrary large or small values of 
io/ek and k/(2kp). We still keep the convergence factor 
exp(— ico0 + ), because that will be crucial in the following, 
where we have to perform the frequency integral over noo 
(~ 1/co for large to). 
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B. Ground— state energy 

The ground-state energy coincides with the grand- 
canonical potential for T — * after adding the term Nfi. 
Using the expression (^), we get (here, np{x) denotes the 
Fermi function at zero temperature) 

E-rpa = ^ n F( e k - A*) e k 

k 

+ oE / — in{i-noo(*,w)[W(fe) + n u (fc s «)]} 



k/0 



2fc^' ~2/7 



) 



(12) 



{Tiuj c — ft 2 j (mbl%)) . Here, the first term is the energy of 
the free CS fermions and the second term is the contribu- 
tion of the magneto-plasmon oscillators. While the first 
sum is confined to the Fermi sphere, one has to discuss 
the convergence of the second sum which we will do be- 
low. In order to simplify the following analysis, we have 
defined 



e(u, v) 




with 



. . 2u 2 f°° dx , M Tl 

e(u, u) = / — m{l + 1 (x, u, v) j 

v Jo n 



(13) 



(14) 



and 



(l-z+) 2 (l-3^) 



(15) 



From now on, we will use A = e 2 fc^/(2/x) as a dimension- 
less measure of the Coulomb interaction. z+(x, u) and 
Z-(x,u) are related to the real and imaginary parts of 
z(x, u) in (|l] 



z±(x,u) 



Ax 2 ± (1 



(16) 



The difference between the energy distributions e and 
e is due to the following: For large to, IIoo(A;,c<;) and 
]lii(&;,a>) decay as w -1 . Therefore, one has to study the 
convergence of the u integral in (|12| ). Here, in the first 
term of an expansion of the ln{- • •} around 1, the factors 
exp(— iujO + ) in Hoo(fc,w) define the integration; in the 
higher terms in the expansion, these convergence factors 
can be put equal to one. e(u, v) is the expression which 
one would obtain setting these factors equal to one ev- 
erywhere. The difference between e and e thus comes 



from the difference in the first term in the expansion of 
the ln{- • •}, calculated once with and once without con- 
vergence factors. 

Now, we analyse the contribution of the magneto- 
plasmon oscillators to -Erpa by studying the dependence 
of e(u,v) on u — k/(2kp)- 

(i) We first discuss the behavior of e(u,v) for large u. 
With 



1 



1 



1 1 - 3a: 



r.2 



u) — x + 0(u~ 2 ) 
a straight forward expansion of e(w, v) yields 



0(u 



(17) 



(I) 2 +2^ 2 



+ 0(u- 5 ). (18) 



Thus, the integral defining -Erpa is convergent for large 
wave vectors k. There is no need of an artificial cut off in 
the ultraviolet. A numerical integration of (M) shows in 
fact that the cut off is rather abrupt, cf. Fig7|j. We shall 
justify below the values v = 1 and A = 0.7, which we 
have chosen here for v and the strength of the Coulomb 
interaction. 



-0.6 



-0.8 



FIG. 1. e(u, 1) for <j> = 2 and A = 0.7. 

(ii) The opposite case of small u leads to a divergence. 
The integral (J14|) is finite; its dominant contributions 



come from . 



\ We find for 



e(0,v) 








(19) 



The result of the numerical integration for general u as 
shown in Fig. ^ confirms a smooth behavior of e(w, v) 
for small u. The negative curvature in the region near 
u = 0.9 is seen to be related to the value of A. This 
feature gets less pronounced as A is decreased. 
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FIG. 2. e(u, 1) for 4> = 2 and A = 0.7. 

Hence, e(u, v) oc — <fi/(8u 2 ) cx k~ 2 for small k and -Erpa 
contains a logarithmic divergence in the infrared regime. 
It is clear from the derivation, that this logarithmic di- 
vergence comes from the first diagram in the RPA series 
in il): 



C. Chemical potential 

The chemical potential, as the ground-state energy, 
should be derived from the grand-canonical potential. 
Calculating the mean particle number N from J^rpa and 
using its value given by the filling factor (||), we get at 
zero temperature 



v = — + 2(j>v 



du u [ud u + 2vd v + Xd\] e(u, v) . (21) 



This equation should determine the chemical potential 
/i via /i = \h~io c lv. But e(u, v) diverges in the infrared 
regime as —<fi/(8u 2 ). Therefore, the r.h.s. of ( ^l| ) be- 
comes infinitely large and there is no positive solution 
for the parameter v from this equation. Hence, we have 
to conclude that within the RPA, there is no satisfactory 
solution for the chemical potential. 



^^X(fc)^n f (k 1 +k)n F (k 1 

k#0 ki 
1~ 2 

(j) [ihuic 



N 



4v 



In- 



i] 



(20) 



Here, we give only the term with the leading singularity 
for decreasing infrared cut off r\\ n <C kp is an upper wave 
vector cut off. The result (^0|) is quickly reproduced if one 
uses only the leading term e(u, v) ~ — <ft/(8u 2 ) in the in- 
tegral for the energy (|l2|). The logarithmic divergence is 
explained as the result of the long-range behavior of the 
CS interaction K(k) in a two-dimensional system. Thus 
we see, that within the RPA, the ground-state energy is 
infinite. 

We wish to add two remarks. First, this divergence in 
the energy and the divergence in the self energy discov- 
ered in Reffl can not have the same origin, since, there, 
the divergence depends on the form of the electron- 
electron interaction, while here, it is already present in 
the absence of a Coulomb interaction and rather a re- 
sult of the pure CS interaction. Second, in other cases of 
RPA calculations, such a divergence in the energy does 
not appear; in the case of the Coulomb interaction in 
three dimensions, for example, one has a convergent in- 
tegral (J d 3 k/k 2 ) and the same is true for the Coulomb 
interaction in two dimensions. Consequently, for the case 
of the CS interaction, it becomes necessary to study dia- 
grams beyond the RPA. 



D. Expansion of the energy in the Coulomb 
interaction 

Here, we wish to study the expansion of the ground- 
state energy in the strength of the Coulomb interaction, 
A. Generally, we have for the energy per particle 



E/N = ~hu c + E^/N + 0(X 2 ) 



(22) 



The first term is the energy of the lowest Landau level; 
in the RPA, it diverges as shown above. Because the 
CS transformation gives an exact reformulation of the 
problem also in the absence of an electron-electron in- 
teraction, this divergence must be compensated by other 
divergent diagrams; we expect s uch a compensation by 
the diagrams discussed in Subsec. Ill E below. In the cal- 
culation of the second term, now, we encounter the prob- 
lem that the unperturbed chemical potential can not be 
determined from the unperturbed grand-canonical po- 
tential because the divergence is not eliminated. As a 
way out, we insert in the second term in ( p2|) the known 
exact unperturbed value /i = huj c /2. Because for free 
composite fermions, \i = \/{2mt,)(AitN / F), i.e. /i = Tiuj c v 
(use eq.(^)) that value of fj, coincides for half-filling ac- 
cidentally with the chemical potential of the free com- 
posite fermions. This substitution fixes the variable 
v, v — 1. Using this and taking the electron densitii 
p w 1.5 • 10 15 to~ 2 from the experiments of Kang et al.El 
in GaAs (e r = 12.8 and m\, = 0.068 m e i), one derives for 
the value of the relative strength of the Coulomb inter- 
action A 0.7 which was used in Figs. 1 and 2 together 
with v = 1. We now obtain for the second term in ( |22] ) 



e 2 20 
l B v 3 / 2 



du 



dx 1 — Z-\-{x, u) 
o 7r 1 +Io(a;,u,t>)_ 



(23) 
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Here, Iq — /|a=o- A numerical evaluation is simple and 
the result is E^p A /N = —0.6 e 2 /ls- This compares rea- 
sonably well with E^/N = -0.466 e 2 /l B obtained in 
numerical simulations in the sphericaLgeometry by Morf 
and d'AmbrumenilEJ, and by Girlichll3 via the threshold 
energy of the many-particle density of states. 

E. Higher order logarithmic divergences 

In the CS problem at d — 2, there are logarithmically 
divergent diagrams in the grand-canonical potential, be- 
cause the CS interaction K(k) diverges as 1/k 2 for small 
wave vectors. In contrast, in the Coulomb problems at 
zero magnetic field in d = 3 and d = 2, such diver- 
gences are absent. The lowest-order divergence is found 
in the logarithmically divergent diagram (|2(]), the first 
of the RPA diagrams. In addition, there are diagrams 
for Q beyond the RPA diverging with higher powers of 
In 77. These diagrams are characterized by independent 
ki -integrations over K{k{); they are either ladder type 
diagrams or maximally crossed diagrams or combinations 
of both. The leading singularities of the ladder diagrams 
in second and third order are 

= ~ ^-Phto^ 4 fin-) , (24) 

= "^ 3) ~ iw^^H) 3 ' (25) 

Similarly, for the maximally crossed diagrams 




The summation of all logarithmically divergent diagrams 
is expected to give a finite result. This should then 
lead to a finite ground state energy and a well defined 
equation for the chemical potential in the absence of an 
electron-electron interaction. Also, such a summation 
should improve the result for the Coulomb contribution 
to the ground-state energy. 



IV. CONCLUSION 

On the basis of bare response functions calculated 
analytically for arbitrary k and uj, we evaluated the 
grand-canonical potential in the RPA without any fur- 
ther approximation. Because of the l/fc 2 -singularity of 
the Chern-Simons interaction and the systems dimension 
d — 2, the grand-canonical potential contains a logarith- 
mically divergent term. Thus, the ground-state energy 
and the chemical potential are not well defined within the 
RPA. The CS transformation leads to the picture of com- 
posite fermions forming a Fermi liquid at zero magnetic 
field, which is experimentally well supported. However, 
for a system of non-interacting electrons in a strong mag- 
netic field, which we call the "unperturbed problem" , this 
transformation already leads to a very complicated refor- 
mulation of the Hamiltonian. We have to conclude that, 
in the CS formulation, one has to go beyond the RPA in 
order to get a satisfying solution of that "unperturbed 
problem" . 

Further, we calculated the first-order term in the ex- 
pansion of the ground-state energy with respect to the 
strength of the Coulomb interaction. There, we used 
for the chemical potential, which remains undefined in 
the RPA, the unperturbed value. The result obtained 
is rs 25% larger than the result of numerical simula- 
tions. In view of the simplicity of the random-phase 
approximation and the absence of a small parameter in 
the "unperturbed problem" , this seems to be a reason- 
able agreement. For an improvement of the calculation 
of the Coulomb contribution to the energy, the "unper- 
turbed problem" must be solved in a better approxima- 
tion, which should be free from a divergence. 

The "unperturbed problem" is well defined, the CS 
transformation is a rigorous reformulation; therefore, one 
should not encounter divergences. Thus, the divergent di- 
agram in the RPA must be compensated by others. We 
finally studied other logarithmically divergent diagrams, 
beyond the RPA, and calculated their leading singulari- 
ties. It is an open question, whether taking these lead- 
ing singularities into account is already sufficient for a 
compensation. We plan to return to this question in a 
following publication. 
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